Zero-point energy of massless scalar fields in the presence 
of soft and semihard boundaries in D dimensions 
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1 Introduction 



In many situations in Quantum Field Theory it is assumed that the fields are defined in a region 
limited by some finite "classical" cavity and submitted to some particular boundary condition. 
"Classical" means here that the boundary which confines the fields has a very precise spatial lo- 
cation and a well defined geometrical shape ("hard" boundary). In the majority of the papers 
discussing the Casimir effect these " hard" boundaries are assumed, although they are unques- 
tionably an idealization. In this scenario, one can argue in what extension does the accumulated 
experience in the Casimir effect depend on the assumption of this kind of boundaries. In some 
papers these conditions have been relaxed — the most recent ones quoted in Refs. [1-5] . 

In this paper we will deepen into the investigation on how the Casimir energy of a massless 
scalar field, defined in a D-dimensional fiat spacetime, depends on the boundary conditions and on 
the dimensionality. For this purpose, three different kinds of confining boundaries are considered: 
"hard" , "soft" and also "semihard" ones. The meaning of this terminology will be clarified latter. 

The problem of determining the expectation value of a physical observable is related to the 
question: how to implement a renormalization scheme in a given situation? In 1948 Casimir 
presented a scheme to obtain a finite result from the divergent zero-point energy of the electro- 
magnetic field 0. Although formally divergent, the difference between the vacuum energy of 
different physical configurations can be finite. If one of these configurations is assumed to have 
a zero vacuum energy, then the difference of the vacuum energy of both configurations is the 
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renormalized one. Therefore the formal definition of the Casimir energy is 

Eren{dn)=Eo{dQ)-Eo{0), (1) 

where Eq^OQ) and Eq{0) are, respectively, the zero-point energies in the presence and in the 
absence of boundaries. In the case of scalar fields, Casimir's approach can be summarized in 
the following steps: a complete set of mode solutions of the Klein-Gordon equation satisfying 
an appropriate boundary condition and the respective eigenfrequencies are found; the divergent 
zero-point energy is regularized by the introduction of an ultraviolet cut-off and, finally, the polar 
part of the regularized energy is removed using a renormalization procedure. 

It is well known that there are two quantities which might be expected to correspond to the 
total renormalized energy of quantum fields [0. The first is called the mode sum energy < E >]^J^'^, 

<^>rn'^=/ rf^-^(iVM-iVoM), (2) 

JO 2 

where is the zero-point energy for each mode, N{uj)duj is the number of modes with frequencies 
between u and u + du in the presence of boundaries and No{uj)dLj is the corresponding quantity 
evaluated in empty space. Eq. (0) gives the renormalized sum of the zero-point energy of each 
mode. The second one is the volume integral of the renormalized energy density, < E >ren^ 
obtained by the Green's function method [Q. In the latter method, in order to calculate the 
renormalized energy for any field, a certain second order differential operator is applied to the 
renormalized Green's function Gren{x,x'), i.e., 

< Tf,y >ren= lim^ -D^z/ [G{x, x') - Go{x, x')] , (3) 



where G{x,x') is the Green's function in the presence of the boundary (dQ) and Go{x,x') is the 
Green's function in the absence of boundaries. Deutsch and Candelas refer to the quantity 
between the brackets as the renormahzed Green's function, since both Green's functions give rise 
to the same uhraviolet singularity structure (as x —>■ x'). If a; belongs to the boundary {dQ) the 
renormahzed stress-tensor < T^,^ >ren can diverge as one gets close to this surface. However, as 
was stressed by these authors, the above argument is not a proof that the renormahzed stress- 
tensor < T^u >ren will diverge as we get closer to (dQ), but it suggests that if the renormahzed 
stress tensor is bounded near (dQ) it means that a delicate cancellation must occur. In the case of 
a perfectly conducting spherical shell in the presence of an electromagnetic field both inside and 
outside the cavity there is a cancellation between the TE and TM modes, giving rise to a finite 
energy density even on the boundary 0. 

It is important to point out that, for the minimally coupled scalar field, such cancellation 
does not occur, which renders the concept of the renormahzed vacuum energy density < Too >ren 
ambiguous. However it is well known that the total renormahzed vacuum energy associated with 
a minimally coupled scalar field obtained by the sum of modes method, < E >™^'^^, must be equal 
to that of the conformally coupled case, since both fields satisfy the same wave equation and have 
the same density of states. Nevertheless, the total renormahzed vacuum energies obtained from 
the Green's function method, < E >^e^, for the minimal and conformal scalar fields, are different. 
Actually, < E >'^°^'' is found to be divergent. Which of these quantities, < E >™°'^'' ox < E >^°^, 
therefore, is the "physical" renormahzed energy of a minimally coupled scalar field? In the bag 
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model this problem is also present [|T0|. Using the Green's function method, Bender and Hays 



11| obtained a quadratic divergence for minimal scalar fields confined in the interior of the bag. 



Also Milton, investigating the zero-point energy of vector fields (gluons) confined in a spherical 
bag |r^|, obtained the same kind of quadratic divergence. 

It has often been suggested by many authors that a full quantum mechanical treatment of 
boundary conditions can solve the above mentioned problem. Recently, Ford and Svaiter have 
confirmed these especulations |^ assuming fluctuating boundaries. By considering confining plates 
as quantum objects with a position probability distribution p, it was shown that this approach 
is able to remove the discrepancy between < E >'^°^^ and < E >^°^ for the minimally coupled 
scalar field, solving a long standing paradox concerning the renormalized energy of the minimally 
and conformally coupled scalar fields. 

There are many other different approaches in order to relax the classical boundary conditions. 
A long time ago, investigating the bag model, some authors discussed quantum corrections to 
this model by quantizing fluctuations around the clasical bag solution [|1^. Working in the same 
direction, Creutz |T^ studied the effects of considering different bag configurations using the 
path integral approach in a theory with different massive scalar fields inside and outside the 
bag. Golestanian and Kardar ^ also use a path integral approach to investigate the problem of 
perfectly reflecting cavities that undergo an arbitrary dynamical deformation. These authors were 
able to calculate the behavior of the mechanical response function (i.e., the ratio between the 
induced force and the deformation field in the linear regime). Some authors, on the other hand. 
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employed a simpler alternative approach, which allows one to deal with more general physical 
situations than the "hard" classical boundary conditions currently used in the literature. They 
imagine a confining "soft" boundary as modeled by a given smoothly increasing potential function 
representing some distribution of matter which interacts with the quantum field |[T|, |^ . Using this 
approach, it is possible to recover "hard" boundaries in some limit. This point will be clarified 
latter. 

The aim of this paper is to discuss the Casimir effect for massless scalar fields subjected not 
only to "hard" boundaries in four dimensional spacetime but also to "soft" and "semihard" ones 
in a general D-dimensional flat spacetime. A classic question to be analysed is what actually 
determines the attractive or repulsive nature of the Casimir force. As it is well known, the sign 
of the Casimir energy may depend on the type of boundary conditions, on the ratio of the flnite 



characteristic lenghts of the cavity and on many other geometrical and topological features |T^, |T6|. 
Recently, the sign of the Casimir energy was discussed in Ref. Jl^ where some results obtained 
in Ref. [0 were generalized, but still assuming only Dirichlet b.c. It is our purpose here to 
address the question of the sign of the Casimir energy in these "new" conflning situations in a 
D-dimensional flat spacetime. 

The paper is organized in the following way. In Section II we review the most simple example 
of the Casimir energy dependence on the ratio between the characteristic lenghts of a "classical" 
cavity. In Section III we analyse the Casimir effect for a minimally coupled scalar fleld in a D- 
dimensional spacetime in the presence of "hard" boundaries. In Section IV we investigate the 
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Casimir effect in the presence of "soft" and "semiliard" boundaries. In Section V we consider tlie 
situation of more tfian one confining potential being imposed onto tlie field. Conclusions are given 
in Section VI. Throughout this paper we use h = c = kB = I- 

2 The Casimir energy in a two-dimensional classical box 

In order to get some insight on the problem of renormalized quantities confined in compact 
domains, we review, in this Section, a well known example. The most simple situation that we 
can imagine in which the vacuum energy is dependent on the ratio of characteristic lengths of a 
cavity is that of a minimally coupled scalar field satisfying classical boundary conditions in a 3D 
spacetime. 

Let us consider a free massless scalar field confined in a 2D rectangular box satisfying Dirichlet 
boundary conditions. Although the presence of the corners are unphysical features, the model 
was also used by Peterson, Hansson and Johnson [|I^ in the study of loop diagrams of a confined 
scalar field in boxes, and it is suitable for our purposes. 

A free real massless scalar field ip{x) defined in a flat spacetime must satisfy the homogeneous 
Klein-Gordon equation. If we restrict the field to the interior of the box, the field modes are 
denumerable and the positive and negative frequency parts form a complete orthonormal set. The 
renormalized energy can be obtained after suitable regularization and renormalization procedure 
of the infinite sum of the zero-point energy of each field mode. Because there is no difference 
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between the density of modes of the minimal and of the conformal scalar fields, the example below 
covers both situations. In this Section we will follow the procedure of Ref . . 

In the Fock representation, there must exist a particular vector |0), called the vacuum or the 
no-particle state. In a 3D spacetime the eigenfrequencies of the field are given by 

'^nma = (^(^-^j +[-J^) J ni,n2 = 1,2,..., (4) 
where Li and L2 are the lengths of the sides of the box. The zero-point energy is 

E{Li,L2) = - J2 ^nin2, (5) 
ni,n2=l 

where ujnin2 is given by Eq. (Q). This expression is divergent and can be written as: 

-i 00 

E^{L,,L2,s) = - E (6) 

for s = — |. 

Eq. (H) is analytic for Re{s) > 1. An analytic regularization method consists in evaluating the 
analytic continuation of the zeta function at the point s = — |. Algebraic manipulations of Eq. 
(I), using Eq. (|), give 

where C{2s) is the Riemann zeta function and A{a, b, s) is the Epstein zeta function defined as: 

00 

A{ai,a2,s)= {ainj + a2nl)~'' . 

n\ ,n2=— 00 



The prime sign in the summation means that the term ni = n2 = is to be excluded. There- 
fore E(^{Li, L2; s) is analytic in the complex s-plane for s G C/{^,1}, and the evaluation of 
Eq{Li,L2, —\) gives the Casimir energy U{Li,L2): 

Instead of analytically regularizing the zero-point energy we can obtain the Casimir energy by 



introducing a suitable cut-off. For the details of these calculations see Ref. |21], and for a general 



discussion about analytic regularization methods used to obtain the renormalized vacuum energy 



of free fields in an arbitrary ultrastatic spacetime see Ref. p2[. A simple inspection of Eq. @ 
shows that the sign of the Casimir energy depends on the ratio between Li and L2, and its behavior 
is shown in Fig. (1) and Fig. (2). In the next two Sections we will extend these calculations to a D- 
dimensional spacetime assuming not only Dirichlet boundary conditions but also other categories 
of boundaries. 
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3 The Casimir energy of a massless scalar field in the 
presence of "hard" boundaries in D dimensions 

Let us consider a free massless scalar field (p{t,x) defined in a. D = d + 1 dimensional Minkowski 
spacetime. If we assume Dirichlet b.c. in a D-1 dimensional box with lenghts Li,L2, . . . ,Ld_i, 
the eigenfrequencies are given by: 

1/2 

(9) 

Using the condition a = Ld_i <^Li,i = l,2,...,D — 2, the energy of the vacuum state is 

1 

(fci)^+...+(w+(^)T- 

(10) 

Note that the summation starts at n = 1 because for the scalar field one should not include the 
modes for which all integers ni, n2, ■ ■ ■ , no-i vanish. As was stressed in the previous Section there 
are two different ways to obtain the Casimir energy using the sum of modes method. The first 
one is to use dimensional regularization in the continuous variable in Eq. (|TU|) and analytically 
extend the Epstein zeta function that will appear after dimensional regularization. A different 
approach is to use dimensional regularization in the continuous variables and to introduce a cut-off 
in the discret one. Let us use this second approach in this case. For "soft" and other types of 
boundaries it may be useful to consider both dimensional and zeta function analytic regularizations 
(see Section IV). 
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ni,n2,...,no-i 



n2'n 



+ . . . + 



The angular part of the integral over the D — 2 dimensional k space can be calculated straight- 
forwardly and if we define the energy per unit area hy = j^i"' ^ — 1; 2, . . . , D — 2, we 



have 



en{a)=F{D)Y: / 
1 Jo 

n=l 



2\ 2 



where 



FiD) 



(2v^) 



2-D 



(11) 



(12) 



The energy per unit area is divergent and should be regularized. Let us introduce in Eq. (|TT]) a 
convergence factor, i.e., an ultraviolet regulator 



exp 



nn 



1 ■ 

2\ 2 



a J 



(13) 



valid for Re{X) > 0. The regularized energy per unit area is finite provided Re{X) > 0, and is 
given by 



oo c 

en{X,a) = F{D)Y: / 



/ nn 
\ a 



1 

2\ 2 



exp 



1 

2\ 2 



(14) 



The Casimir energy per unit area (the renormalized vacuum energy per unit area) is defined by 



UD{a)= lim [eDiX,a) + eD{X,R-a)-eDiX,r]R)-eDiX,{l-7])R)], (15) 

A— >0,iJ^co 



where rj is a real number between zero and one. A straighforward calculation gives (see ||2l[] ): 

1 r(f)c(D) 

UD{a) = --—IT — — • (16) 



(47r)f a^'^ 
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Thus the Casimir energy is negative for any D in this particular configuration. This result is in 



agreement with Ambjorn and Wolfram |jT6|. After a schematic review of this well known case, we 
are now in position to investigate two different kinds of boundaries: the "soft" and "semihard" 
ones. It is important to stress that since we are using the sum of modes method to find the 
Casimir energy, both the cases of minimally and conformally coupled scalar fields are covered. As 
we discussed before, this comes from the fact that there is no difference between the density of 
modes of the minimal and of conformal scalar fields. 

4 The effect of "soft" and "semihard" boundaries in the 
Casimir energy 



In this Section we will investigate the Casimir effect of a massless scalar field in the presence of 
"soft" and "semihard" boundaries in a general D-dimensional spacetime. The idea is to replace the 
"hard" Dirichlet walls by some confining potential V{x) in the xd-i direction (the Dirichlet 
condition corresponds to the particular case where V{x) vanishes inside the cavity and V{x) 
becomes infinite on the boundary). In this case, the spatial modes of the scalar quantum field 
satisfy a Schrodinger-like equation and its spectrum will be denoted by u^. This confining potential 
may be interpreted as representing some distribution of matter with which the quantum field 
interacts. Since the potential acts as effective plates, it is relevant to state that the Casimir forces 
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will act upon the matter distribution modeled by V{x). When all modes are completely supressed 
only for x oo, the effective boundary is called "soft". We can also imagine an intermediate 
situation between "hard" and "soft" — that may be called "semihard" — where the complete 
supression happens for a given finite x value. In this case, the potential V{x) decreases smoothly 
from an infinite value on the boundary surface dQ to = far from dQ. In this sense, Actor and 
Bender atribute a sort of "texture" to the boundary effective surface [|[], where the case of the 
harmonic oscillator potential in a particular direction, say xd-i, was investigated ioi D = 4 (see 
also i). 

Assuming that the boundary conditions in the direction xd-i = x are dictated by a generic 
potential V{x,a), (where a is a characteristic length of the system), the vacuum energy per unit 
area can be written as 

1 ^ fl-D\^, D-l 



2(47r)— V 2 / „ 

The first situation that we would like to discuss is that of a potential which is "semihard" 
near the origin and "soft" for large x. An example of such situation can be given by the following 
potential, plotted in Fig. (3): 



where a and h have dimension of \length] and Vq has the dimension of \length] ^. The solution of 
the Schrodinger's equation in the case 6 = a is well known ||2^. However, the above potential is 



more suitable to work out some limits and indeed only slight changes are needed to get the solution 
for h ^ a. It is straightforward to show that the energy levels of the Schrodinger's equation for 
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h ^ a are given by 



V 62 



1 1 
n H h - 

2 4 



n = 0, 1,2, 



(19) 



Substituting Eq. (|T9[) in the vacuum energy density given by Eq. (|T7[) we obtain: 



ez)(a, 6) 



1 



V 



2(47r)V V 62 y - V 2 
After doing the analytic extension of the Hurwitz zeta-function 



c 



(20) 



C(^;g) = E 



1 



{n + qY' 

which is analytic at the beginning of an open connected set of points of the complex plane, i.e. 
Re{z) > 1, we obtain for the vacuum energy density: 



eD{a,b) 



1 



( 



D-l 

32Vn\~ 



e*-l ^0 n\ n- 2 



2(47r)V V 62 

In the above equation Bn{q) are the Bernoulli polynomials p4[ and q is given by 



D+l 



(21) 



(22) 



To obtain the Casimir energy we have to subtract the polar part of the above equation, which 
is easily seen to be a single term in the summation, since the integral is finite. Two comments 
are in order: the first is that a renormalization procedure is necessary only for odd- dimensional 
{D = 2m — 1, m = 1, 2, . . .) spacetimes, because Eq. (^) is already finite for even D; second is 
that the "soft" boundaries change the structure of the poles of the model, i.e., the residues of the 



13 



polar part are given by {D is an odd-integer number): 

Res {eD{a,b)) = Bp+i^q). 

Two limits are of interest (we keep h fixed from now on): (i) a? <^ (Vo)~^ (hence q 3/4) in 
which case the potential behaves just as a "hard" impenetrable (Dirichlet) wall for x ^ 0, while 
it behaves as a harmonic oscillator potential restricted to x > 0, and (ii) ^ (K))^^ (? V^)- 
In the formalism of Actor and Bender it is also possible to obtain limit (i) from the harmonic 
oscillator (HO) result by just discarding some of the eigenvalues in the zeta-function, and this can 
be called the \H0 limit. 

First let us investigate the limit ^ (Vo)~^. In order to compare the Casimir energy with 
the value obtained in Ref. [|1|, where the = 4 case was treated, we need to use the particular 
value of the Hurwitz zeta-function C(— 3/2; 3/4) = 0.02093. In addition, one can still define 
A = (Vo/&^)^^^^, with dimension of [length], and interpret A as the "characteristic separation 
distance" between the "hard wall" at x = and the "soft" one at a; ~ A. In this way, from Eq. 



(|20|) we readily obtain: 



27/4 

e4(a, b) = - A-3C(-3/2; 3/4), (23) 

and £4 is half the value found in Ref. It should be noted, however, that the value found in 
Ref. for the well known Casimir energy in D = 4 between two Dirichlet plates, as the limit of 
their result for the harmonic oscillator potential, is also twice the value found in Refs. |IB|, [TB|. In 



the limit of large separation between the "walls" (A 00), one recovers the free half-space result: 
e 0^. 
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Let us examine now the values of the vacuum energy density given by Eq. (pTD, for D = 2, 3, 4, 
in both hmiting cases ^ and <^ (K))^^ (o' 1/2 and g — 3/4 respectively). An 

interesting feature for low dimensional spacetimes, i.e., for D = 2 and D = 3 is found: indeed, 
when ^ the vacuum fluctuations give rise to a repulsive force corresponding to energy 

densities 

^lim 62(a, 6) = [+0.0724]A-^ (24) 

and 

lim €3(0,6) = [+0.0160]A~^ (25) 

respectively, while in the limit <^ (M))^^ the force becomes attractive, corresponding to energy 
densities 

hm e2(a,6) = [-0.0562]A-i (26) 

and 

hm e3(a,6) = [-0.0108]A^2_ ^27) 

Thus, for D = 2,3, there must exist some finite a for which the Casimir energy vanishes. The 
same behavior was found in the case of the two-dimensional classical box considered in Section II, 
where the sign of the Casimir energy was shown to depend on the ratio between the lengths of 



the box |17|. It is important to emphasize that single poles appear in both limits only for D = 3. 
In the four-dimensional case, for both limits, the Casimir force is found to be always attractive. 
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with corresponding values for energy densities 



lim £4(0,6) = [-0.0059]A"^ (2^ 

a2>(Vb)-l 



and 



lim e4(a,6) = [-0.0075]A^^ (29) 

and the latter is exactly the value of Eq. (pS]). Only in the limit a? <^ (Vo)~^ the Casimir energy 
has the same (negative) sign as in the lowest spacetime dimensions, while in the limit ^ (K))~^ 
the sign of the Casimir energy for D = 4 is opposite to those oi D = 2 and D = 3. 

In Table 1 we present the values of the Casimir energy for a massless scalar field in the presence 
of the asymmetric potential V{x,a,h) in both limits a? <^ (Vq)^^ and a? ^ (Vq)"^, for spacetime 
dimensionality varying in the range 2 < D < 20. Although we went up to D = 100, we report 
on the table, for simplicity, only the values for D < 20. From these numerical results we can 
discuss the sign of the Casimir energy, which is not straightforward from the analytical expression 
Eq. (^). The first point we would like to stress is that, independently of a, the Casimir energy 
€£) has the sign (— l)''"+i for D = 4(m + l),m = 0, 1,2,3. For m > 4, the sign becomes (— l)"*; 
indeed, for D > 18 this pattern of signs is broken. In any case, we have checked that even for the 
wide range 2 < D < 100 the modulus of always decrease (we have assumed always A = 1). 
Secondly, we would like to point out that, for this potential, it is possible to conceive a gedanken 
experiment to investigate the possibility of higher dimensions. Whenever it is found, varying a, 
that the Casimir energy vanishes (and hence changes sign), it is safe to assert that D = 41 + 2 
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OT D = Al + 3 {I being a natural number). Alternatively, as for D = 4, if the sign never changes, 
spacetime dimensionality should he D = Al or D = 41 + 1. This is the maximum we can conclude 
from the experience in this case. But once again we stress that the sign pattern is altered for 
D > 18: while for D < 18 and e^+i have the same sign, for D > 18 they have opposite signs. 

The second situation that we would like to discuss is the case of an increasing potential in the 
xd~i direction which becomes infinite for x = and x = a. To represent this situation let us 
assume that the potential is given by 

V{x) = V^'^ cot^ {^x^ . (30) 

Using the Actor and Bender terminology [|ll, one can say that this situation is equivalent to two 
"semihard walls", one at x = 0, the other oX x = a. Note that, in this case, the field is confined 
to the interior of the region < a; < a (see Fig. (4)). It is straightforward to show that the values 
of z/^ are given by 

i^l = %{n' + An(3-2(3), n = l,2,... (31) 
a 



where j3 is defined by 



It follows that the vacuum energy density is now given by 

eoia) = -^^r (^-j + 4n(3 - 2/?)^. (33) 

Here it is necessary to analytically extend a modified Epstein zeta-function. This was done by 



Ford and also by Birrell and Ford . In this case the analytic expression for the vacuum energy 
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density is very complicated and it is not reported here, but an important difference between the 
first potential and this second one is that in some limits it is possible to recover exactly the 
Dirichlet walls, i.e., perfectly conducting plates separated by a distance a. These limits can be 
obtained by expanding the potential around the point x = ^ and neglecting terms higher than 
the second order. Let us first consider the limiting case Vq — 0. We see that, in this case, we 
come back to the problem of classical parallel plates (Dirichlet b.c.) placed in the xd-i direction. 
When Vq ^ also (3^0 and the Casimir energy reduces again to (see Eq. (0)): 



l_r(i)CP) 

(47r)f 



^D{a) = ^ (34) 



after making use of the reflection formula for (^(1 — D) 



T{D/2)n-^/\{D) = r C--^) vr^C(l " D). 



For large values of the eigenvalue /?, and for n <ti P, i.e., for the lower levels, we get: 

ul = ^in + l/2). (35) 

This limit is completely analogous to the case of the harmonic oscilator potential studied in Refs. 
[]1], 0. In the next Section we compute the vacuum energy of a scalar field in some configurations 
where the field is constrained by different potentials in different directions. 
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5 The Casimir energy of a massless scalar field in a hy- 
perbox with diff^erent boundary conditions 



In Section IV the scalar field was supposed to be constrained by a hyperbox where only in one 
direction the "hard" Dirichlet plates were replaced by a confining potential. In this Section we 
analyse different situations in which, out of the D dimensions of spacetime, in p of them the field 
should satisfy Dirichlet boundary conditions, and in each of the remainder D — p — 1 directions 
it is subjected to confining potentials. Moreover, as we are working in Cartesian coordinates, 
our formalism allows us to choose different potentials acting upon the field in each one of the 
remainder D — p — 1 directions. 

In the directions x*, i = 1, 2, ..,p, we impose the vanishing of the field at parallel plates located 
at — 0,Li,i — 1,2, ... ,p (Dirichlet boundary conditions); in the directions ,j — p + l,p + 
2, . . . , D — 1, we choose D —p — 1 potentials Vj{xj, aj) (all of them may be a priori different), each 
one depending upon different characteristic sizes Uj. The vacuum energy is easily written as: 



Now, if all Li are made much greater than all aj, then we can replace the first p summations 




1/2 



(36) 



where the functions 



.{ttj) represent the spectra of eigenvalues of the corresponding Schrodinger's 



equation. 
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by integrals: 



dPk 



j=p+i 



1/2 



(37) 



^ \i=l / np+i,...,nD-i (27r)^' 

where k is the length of the vector {ki, .., kp). The integral above is in a well-suited form to apply 
dimensional regularization, with the result: 



E 



D 



2(47r) 2 Vj^i y V ^ / np+i,...,nD_i 



;aj,+i) + ... + i/^^_^(ai)_i) 



i+p 



(38) 



Instead of discussing the general case, let us work out two specific cases, considering a 4D 
spacetime. In one spatial direction (for example z) let us impose Dirichlet boundary conditions 
(with plates separated by a distance L); hence p = 1. Besides, we will subject the field to the 
same potential in the x and y directions: 



^(..a.6)^^o-(My 



(39) 



and 



y(y,a,6) = yo 



1/2 / a y 



The respective spectra of the Schrodinger's equation are already known to us: 



2 f^^VoY 



(40) 



(41) 



and 



/32^ 
V 62 



(m + q), 



(42) 
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where the value of q is given by Eq. (0). Exploiting the symmetry between x and y directions, 
the summation which then appears from Eq. (B^), 



^ [n + m + 2g] , 

n,m=0 

can be put in a more tractable form by using (see p6|): 



(43) 



^ [n + m + c]-' = C{s - 1; c) + (1 - c)C(s; c). 

n,m=0 



(44) 



From the general result of Eq. (pHD, we worked out the limits a? ^ (Vq) ^ and ^ (Vq) ^ 
(respectively g — > 1/2 and g ^ 3/4). In the limit a? ^ (Vo)~^ the vacuum energy is given by: 



a2>{Vb)- 



dt4 — + y 



(45) 



It is now easy to obtain the energy density, i.e., energy per unit area, from the expression above. 
First one divides the expression by L. Although the plates in the x and y directions are replaced 
by confining potentials, one can still assign to these directions a "characteristic distance" between 
the "walls", as stressed in the previous section, given by A = ^ . Therefore, in order to 

obtain the energy density e, one further divides Eq. (|45| ) by A, which yields: 

1 



_ ^ 

lim eAa, b) = \ ^ — 

a2>(yo)-l TT 



oo f-3 oo ID 



(46) 



where the integral is finite and the polar part, given by the fourth term in the summation, is 
identically zero because -B3 = (otherwise it would be discarded as usual). In this way, the 
Casimir energy in this limit reads: 



lim €4(0,6) 
a2»(yo)-i 



-0.0069]A 



-3 



(47) 
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which is quite the same value of the hmit obtained in Section IV, Eq. ( pH]) (confining potential 
only in one direction). 

In the other limit <^ (K))^^ the vacuum energy density is given by (see Eq. 

lim eJa,b) = A ^ — 

a2<(Vb)-l TT 



1 e^-l ^r'n n\ n-3 ^ 



n=0 

+ 1 I r dx '"'^'''^'^ I V (-l)"^-B>n(3/2) 1 ^ 

where each summation contains a pole, with corresponding residues ~^'-^'^/'^'> = ^ and ^^^^^'^^ = 11. 
The regularized vacuum energy density or, simply, the Casimir energy in this case is evaluated to 
give 

lim €4(0,6) = [-0.0310]A~^ (49) 

also negative, giving rise to an attractive Casimir force. This value is to be compared to that 
of Eq. (pQl). Thus, in any case, the replacement of two parallel Dirichlet plates in one further 
direction, y, in comparison with the case of Section IV, makes the absolute value of the Casimir 
energy to increase. 

For completeness, let us calculate the Casimir energy when this "soft" potential acts in all 
three spatial directions, again for D = 4; so p = in this situation. We can obtain, in this case, 
the vacuum energy density from Eq. (^), by simply dividing it by = (^^j ; it reads: 

e^{a,b) = +2'/^\-^ J2 [n + m + l + Sq]'/^. (50) 

n,m,l=0 
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Again this summation can be simplified EH] if use is made of the relation: 



°° 13 1 

E [n + m + l + c]-' = -a-s - 2; c) + (- - c)C(s - 1; c) + -(c - l)(c - 2)C(s; c). (51) 

n,m,(=0 

In the limit ^ (^)~^ the Casimir energy is given by: 

lim eJa,b) = \-0.0l32]\-^, (52) 
and in the limit <^ (K))~^ it reads: 



lim €4(0,6) = [+0.0255]A"^ (53) 

For this new configuration the result is qualitatively different from the situation with confining 
potentials only in two directions (previous case). In the limit ^ (Vq)"^, €4 is still negative and 
twice the value of Eq. (^7D; in the other limit, <^ (Vq)"^, €4 changes sign but its absolute value 
is smaller than the value of Eq. (^9]) . 

6 Final remarks 

We have examined how the Casimir energy of a massless scalar field confined to the interior of 
a D — 1 dimensional hyperbox depends on different kinds of boundary conditions and on the 
dimensionality of spacetime. "Classical" Dirichlet boundary conditions in one, two and three 
directions were relaxed; in these directions the constraints on the field were assumed to be given 
by some smoothly increasing potential that represents some distribution of matter which interacts 
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with the quantum field. The new contribution of the present work regards the study of the Casimir 
effect generated by two different types of boundary conditions, namely, the "soft" and "semihard" 
ones. In particular, we have discussed in details the case of a confining asymmetric potential 
V{x) = VQ^'^{a/x — x/by, which presents the feature of being "soft" for a; oo and "semihard" 
for X ^ 0. Although the choice of the potential is in general dictated by the solvability of the 
Schrodinger's equation and by the manageable structure of the sum of proper modes, the study of 
how the Casimir effect depends on the boundary conditions opens new perspectives, which could 
lead to a deeper understanding of the interaction of real (not perfectly conducting) boundaries 
with the field. 

Let us stress now some remarkable differences between the case studied here and the case of 



the Casimir energy in a D-dimensional "hard" hyperbox considered in Ref. |18|. The first one 



regards the repulsive or attractive nature of the Casimir force. In Ref. [0] it was shown that 
the force is attractive if the number p of finite and equal edges of a rectangular box is odd or for 
very large even values of p, irrespective of D. However, it was also shown in that paper that for 
each small evenp there exists a critical spacetime dimension Dc{p) such that the force is repulsive 
if D < Dc{p) and attractive otherwise. Our calculations have shown that, for the asymmetric 
potential considered here, there is no critical dimension Dc such that for D > Dc the Casimir 
energy have always the same sign. What we have demonstrated is that, independently of the 
parameter a, there is a regular pattern for the sign of the Casimir energy for D < 18; for other 
values of D this sign pattern is broken (see Table 1). The computation of the vacuum energy up 
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to D = 100 showed that {eol always decrease with increasing D. 

A second comment concerning the sign of the Casimir energy is related to how it changes when 
p' pairs of "hard" Dirichlet walls are substituted by "soft" and/or "semihard" potentials. It is 
well known that the "hard" wall Casimir energies in D = A are negative for p = 1,3 and positive 
for p = 2. A significant difference between the "hard" and the "soft" wall Casimir effects for long 
waveguides was first pointed out in Ref. and occurs only in the case where two pairs of Dirichlet 
plates were replaced with harmonic oscillator potential; when one or three pairs were substituted 
by the same potential, there is a qualitative similarity between the two cases. Our result, based 
on a different smoothly increasing potential, corroborates this tendency, but attention should be 
drawn to the following point when p' = 3. Our result in the limit <^ (^o)~^ should be compared 
to what is called the ^HO limit of Ref. which is positive and equal to our result (up to a 
systematic factor of 2, as stressed in the text). However, in the other asymptotic limit ^ (Vo)~^, 
the sign changes. Similarly, regarding "hard" walls for D = 4 and p = 3, it is known that in the 
symmetric configuration with equal edges (Li = L2 = L3), the Casimir energy is negative, while 
in Ref. |T^ it was demonstrated that allowing Li 7^ L2 7^ L3 the sign can also change. 



As a last comment let us revisit the question of the instability of the semiclassical Abraham- 
Lorentz- Casimir model for the electron. The original Casimir's idea was that the electrostatic 
repulsion due to the electron's distribution of charge could be balanced by the zero-point fluctu- 
ations of the electromagnetic field inside and outside the conducting shell, assumed at that time 



to be "hard". Unfortunately, in D = A, the Casimir force is found to be repulsive [27]. How this 
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fact depends on D was analysed in Ref. [|T8]. There, it was shown that the stabihty of a Casimir 
electron model in higher-dimensional spacetimes would be possible only for a number of dimen- 
sions D > 8. Therefore, one can imagine a toy model of a stable semiclassical electron where the 
Poincare stress has quantum electromagnetic origin only if one lived in a higher- dimensional flat 
spacetime. All these results take into account that the boundaries of the electron are "hard" ones. 



In Ref. IjT^ it was shown that a negative zero-point energy can be obtained for such b.c. only 
for a very unexpected particular (and antisymmetric) shape and size. Nevertheless, in the light of 
the new results for confining "soft" boundaries, the Casimir idea of how the semiclassical electron 
could be stabilized may be revived. Indeed, there are two results that could be interpreted as an 
indication in this direction. First, it was demonstrated in Ref. ^ that the Casimir energy of a 
spherical "soft" cavity is negative. Second, Eqs. (|52D and (|53|) show that it is possible to find a 
particular value of a which compensates the electrostatic repulsion in D = 4. Thus, these results 
suggest how the hypothesis of a perfect conducting shell confining the electron was overwhelming. 
On the other hand, these results give rise to an important general question which, to the best of 
our knowledge, has no general answer yet, namely, how the sign of the Casimir energy changes 
when one changes the physical parameters and boundary conditions. It seems to us that the only 
way to discuss the attractive or repulsive character of the Casimir effect is, up to now, by direct 
computation case by case. 

Finally, it should be pointed out that the mode sum energy defined by Deutsch and Candelas 
stresses the fact that a classical "perfect conductor" boundary condition is unphysical and there 
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is a sufficiently high frequency uq for which the modes are not confined by the plates (in the case of 
dieletric materials this is called the plasma frequency). In view of this, the modes in the continuum 
will cancel out, leading us to assert that the only relevant modes to consider in the Casimir effect 
are the discrete ones. A natural extension of this paper is to consider a partially transparent 
boundary, which can be modeled, for example, by the modified Poschl- Teller potential given by 
V{x) = ~^ce cosh^^ax - Another direction to look on is to investigate these confining potentials in 
different geometries, as for example a spherical one, trying to generalize the results obtained by 



Bender and Milton p8| . Interesting physical situations are those of a partially transparent sphere 
and spheres with "soft" and "semihard" boundaries. It is clear that this problem is of great interest 
in the framework of the bag model and may shed additional light on the Abraham-Lorentz-Casimir 
model for the electron. This subject is under investigation by the authors. 
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Table I: Some limits of the Casimir energy for confining potential a, 6) = Vq (f ~ f j 
one direction and Dirichlet b.c. in the other D — 2 directions for various D, where A = (V^)/6^)-i/^. 




Figure 1: The Casimir energy of a two-dimensional classical box U{Li,L2) as a function of its lengths. 
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Figure 3: The asymmetric potential V{x,a) = V^^"^ (| — |)^, which is "semihard" for x ^ and "soft" 
for X — > oo. 
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Figure 4: The semihard potential V{x) = V^^^^cot^ (f^)' which reduces to Dirichlet walls separated by 
in the limit Vq —>■ 0. 
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